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Abstract 

We study the effects of dissipation and decoherence induced on a harmonic oscillator by the 
coupling to a chaotic system with two degrees of freedom. Using the Feynman- Vernon approach 
and treating the chaotic system semiclassically we show that the effects of the low dimensional 
chaotic environment are in many ways similar to those produced by thermal baths. The classical 
correlation and response functions play important roles in both classical and quantum formulations. 
Our results are qualitatively similar to the high temperature regime of the Caldeira-Leggett model. 
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I. INTRODUCTION. 



The relation between chaos and the phenomena of quantum dissipation and decoherence 



I quantum aissipation 



IG 
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The 



has attracted a lot of attention in the last ten years 
problem considered in most works involves the weak interaction of a chaotic system with 
an external oscillator. Various points of view have been considered by different authors. 
One approach concentrates on the chaotic system itself, focusing on the dissipation and 
treating the external oscillator as a time dependent parameter that perturbs the chaotic 
system Q, Q, |^ . The basic assumption is that the external system is slow and sufficiently 
heavy not to be affected significantly by the coupling 0, 0, 0, Q . The effect on the chaotic 
system, on the other hand, is that of an adiabatic perturbation. Under these conditions, a 
dissipative force, acting on the external system, may result. 

A different view of the same problem focus on the semiclassical limit of chaotic systems. 
It Ms been shown gfly flfi that the coupling of a chaot.c .y.te,n with an external 
environment, represented implicitly by a small diffusion constant in the classical and quan- 
tum versions of the Focker-Planck equation, leads to a very close correspondence between 
the classical and quantum evolutions. The coupling causes destruction of quantum interfer- 
ence and, at the same time, it washes out the fine structures of the classical distributions, 
bringing the two dynamics together. 

In this paper we consider the chaotic system and the external oscillator explicitly, as a 
single globally conservative Hamiltonian system. We use the Feynman- Vernon approach to 
trace out the chaotic system variables and construct an effective dynamics for the oscillator. 



in close ana! 



Leggett 



ogy with the treatment of the Brownian motion considered by Caldeira and 



14| . While we focus on the oscillator, looking at dissipation and decoherence, 
the effects of the coupling on the chaotic system are also taken into account consistently. 
This characteristic of the Feynman- Vernon method turns out to be very important in this 
problem, since both the oscillator and the chaotic system have small number of degrees of 
freedom and are both affected by the mutual interaction. However, whereas the treatment 
of Caldeira and Leggett is phenomenological, in the sense that the spectral properties of the 
reservoir are not derived from its Hamiltonian, the case of a small chaotic environment has 
to be treated dynamically. In other words, dissipation and decoherence have to come out 
directly from correlations and response functions. 
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Our purpose here is to understand under what conditions a chaotic system with only two 
degrees of freedom can produce dissipation and decoherence, phenomena usually related to 



many body thermal baths. In a previous paper jl5| we have considered the interaction of an 
oscillator with a chaotic system from a classical point of view. We showed that the effects 
of the oscillator on the environment cannot be neglected. Here we consider the quantum 
version of the same problem, assuming the chaotic system to be in the semiclassical regime. 

In treatment of the Brownian motion by Caldeira and Leggett, the degrees of freedom 
playing the role of the environment are averaged with a canonical ensemble, since the reser- 
voir is kept at the constant temperature. Here, since the environment is small, the micro- 
canonical ensemble is more adequate. A similar approach was recently considere d by Espos- 
ito and Gaspard using random matrix theory to model the chaotic environment |l6l . 

The paper is organized as follows: in Sec. II, we review some aspects of the classical 
formulation that are useful for the quantum analysis. In Sec. Ill, we present the quantum 
formulation using Feynman- Vernon approach The formal development leads to quan- 
tum correlation and response functions, that we calculate semiclassically. In Sec. IV, we 
consider two basic applications: first, the propagation of an Gaussian state, where we char- 
acterize quantum dissipation. Second, we calculate the evolution of a superposition of two 
Gaussian states, focusing on the decoherence due to the chaotic environment. In section V 
we present our conclusions. 

II. CLASSICAL FORMULATION 

In this section we describe the behavior of a system of interest coupled to a small chaotic 
environment from the classical point of view. Although the formalism presented here can 
be extended to more general systems, we particularize right away to the case of a harmonic 
oscillator interacting with the so called Nelson system. The discussion outlined here is a 



summary of the detailed results presented in ref. 
of the system is given by 



(see also 



mm 



). The Hamiltonian 



H = H,{z)+H,{x,y) + Vj{x,z), (1) 

where 

^.W^l,.^. (2) 
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represents the system of interest, 



V/(x, z) = •yxz, 



(3) 



is the interaction potential and 



Hc{x, y) 




(4) 



is the chaotic Hamiltonian, known as Nelson system (NS) j20|. The NS exhibits soft chaos 
and is fairly regular for £'c^0.05, strongly chaotic for £'^^^0.3 and mixed for intermediate 
values of the energy. 

In order to investigate the situation where He plays the role of an external environment for 
the oscillator, we assume that detailed information about the chaotic system is not available. 
If the environment were modeled by a heat bath, the only macroscopic relevant information 
would be its temperature. In the present case we assume that the only information available 
is the initial energy -Ec(O) of the chaotic system. For the oscillator this implies that only 
averages of its observables (over the chaotic system variables) are accessible. 

When the coupling between the chaotic system and the oscillator is turned on, the overall 
conserved energy flows from one system to the other. The oscillator's energy, in particular, 
fluctuates as a function of the time for each specific trajectory. The oscillator's average 
energy is calculated by taking an ensemble of initial conditions uniformly distributed over 
the chaotic energy surface -Ec(O). For the oscillator we fix only one initial condition, which 
we choose to be zifd) = and Pz{0) = >/2mEjfi). The microcanonical ensemble of chaotic 
initial conditions plus the fixed oscillator's initial condition are propagated numerically and, 
at each instant, Hn is calculated for each trajectory and its average value is computed. We 
have shown in [15[ that the oscillator's average energy {Eo{t)) tends to a constant value for 
long times, indicating a 'thermalization' of the systems. 

The short time behavior of {Eo{t)) can be obtained from the Linear Response Theory 
[2^ . From the equations of motion for z and Pz we find 





(5) 



(6) 



where Zd(t) and Pzai^) are the decoupled solutions, given by, 

pJO) 

Zd{t) = — — sin (uot) pz^it) = pz{0) 




(7) 
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where r(t — s) = sin [u!o{t — s)]/u!o and x(t — s) = cos [ujoit — s)]. Thus, 

^0 

7^/ ds / dux{t - s)x{t - u){x{s)x{u)) (8) 
7o ^0 

and 

(z'it)) = zlit) - '^-Izait) 1^ dsTit - s){xis)) + 

^ [ ds [ dur{t- s)r{t-u){x{s)x{u)). (9) 
^ Jo Jo 

The oscillator's average energy can now be obtained from 

Equations (jH)) and Q show that we need {x(t)) and (x(O)x(t)) in order to calculate 
{Eo{t)). The calculation of such averages involve the distribution function p{q,p]t) whose 
initial value is p{q,p; 0) = 5(i/c(g,p)-^c(0))/S(E,(0)), with S(E,(0)) = / dqdp6{Hc{q,p)- 
Ec{0)). Here we are using q = (x, y) and p = {px,Py) for the coordinates and momenta of the 
chaotic system. If the chaotic system were isolated, p would be an invariant distribution and 
p{q,p;t) = p{q,p;0). The coupling, however, causes the value of Hc{q(t) , p(t)) to fluctuate 
in time, distorting the energy surface = Ec{0). Linear Response Theory provides the 



2l|. Keeping in 



first order corrections to this distribution in the limit of weak coupling 
and Q only terms up to order 7^, we find 

(x(t)) = {x{t))e - 7 / ds(f)xx{t - s)z{s) (11) 
Jo 

and 

{x{0)x{t)) = (x(0)x(t))e, (12) 

where {A{q,p))e = f dqdp A{q,p)pe{q,p) and pe = S{Hc — Ec{0)) /T,{Ec{0)) is the micro- 
canonical distribution of the isolated chaotic system, (pxxit) is the response function, given 
by H 



0...(t) = ({a;(O),a;(t)})e= / / d\/pe {x(0), (13) 
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where dV = dx{0)dy{0)dpx{0)dpy{0) and {.,.} is the Poisson bracket with respect to the 
initial conditions. Since Hc{—x) = Hc{x), {x{t))e = 0. Substituting (fTT| and (fT^ into ^ 
and we obtain 



(Plit)) = plit) + 27V.,(t) /o dsxit - s) /; dM0,.,.(s - u)zd{u) 
+7^ /o '^■^ /o '^^^(^ ~ ■^)^(^ ~ n)(x(s)a;(M))e, 



(14) 



(15) 



{z\t)) = zj{t) + ^^z,{t) dsT{t - .) /; du^Us - u)z,{u) 

+ i lo lo ^"^(t - s)T{t - U){x{s)x{u))e. 

Eqs. (fT^ and show that all the influence of the chaotic system is contained in the 
functions (x(0)x(t))e and (pxxit). For NS, the response function is given by j3| 

0xx(t) = -^(p.(O)x(t))e. (16) 

The correlation functions {px{0)x{t))e and (x(0)x(t))p are obtained numerically with a flxed 
time step symplectic integration algorithm [22| applied to the isolated chaotic system. Fig. 
1 shows the numerical correlation functions for Ec{0) = 0.38. These numerical results can 
be well fitted by the expressions 

{x{0)x{t))e = ere""* COS cut, ^^^^ 

with decay rates a = 0.0418 and (3 = 0.0456, amplitudes a = 1.865 and fi = 0.409 and 
frequencies of oscillation uj = 0.1963 and Q = 0.2043 with ~ 10~^ (see Fig.l). Notice 
that the exponents a and f3 and the frequencies u and Q are very similar. If the interacting 
system were integrable, these functions would exhibit quasi-periodic oscillations. 

Considering the expressions ()17|) and assuming Q ^ u and (3 ^a,we obtain the following 
result for (^o(i)) 

{EM = EM + 2!(S + At + fit) + git)), (18) 
m 

where i? is a constant, fit) is an oscillatory function and git) is proportional to e~°*. The 
important result is the coefficient A 



^i.l + .- + a-)-^^ 



[iuj^ - uj)^ + a-'Wiuj^ + u)^ + 



(19) 
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For fixed oscillator frequency ujq and a given chaotic energy shell Ec{0) (and, consequently, 
for given a, fi, u and a), the ratio Eo{0) / Ec{0) is the responsible for the average increase or 
decrease of {Eo{t)). The short time equilibrium in the energy flow is given by the condition 
A = 0, or 

-{u^ + u^ + a^). (20) 



Ec{0) AfiLU 

The equation of motion of z{t) under the average effect of the interaction with the chaotic 
system can also be written in terms of the response function as 

m) + ^oim) = --m) = - f dsct^ut-s){z{s)). (21) 

m m In 



Integrating by parts yields 

m) 

where 



Ll-I^IM] {z{t)) + ^- fdsFit-s){zis)) + ^ziO)F{t) = 0, (22) 
\ m J Jo ^ 



f 2/ie""'^*"*^ 
F{t -s) = ds4>^^{t -s) = — ^{wcos [uj{t - s)] + asin [uj{t - s)]}. (23) 



Eq.()22|) shows that the interaction produces a harmonic correction to the original poten- 
tial, a dissipative term with memory and an external force proportional to z{0). The choice 
-2(0) = simplifies (j^^ and turns it into an average Langevin equation. 

Fig. 2 shows a comparison between the numerically calculated 'bare' oscillator energy 
{Eo{t)), where 

^oW^i^^, (24) 

the 're-normalized' oscillator energy {Eor{t)), where 

and the expression (jlSj) without the oscillating term f{t). We have chosen 7 and m so that 
uJq — 7^F(0)/m > 0. We also have chosen cuq so that e~"/'^° ~ 10~^ and g{t) decreases very 
fast. In this case only the linear and the oscillating terms in Eq. ()18|) are important. We have 
subtracted the oscillating part of Eq. ()18j) in Fig. 2 to highlight the linear increase or decrease 
in the average energy. In the time scale of Fig. 2, which corresponds to several periods of the 
decoupled oscillator, the linear behavior describes very well the numerical results. Fig. 2(b) 
shows the equilibrium situation according to Eq.()20p. Notice that F{t) decays very fast in 
the time scale I/cuq, leading to the dissipative force in Eq.( 



In the next sections we consider the quantum counterpart of these classical calculations. 
The chaotic system will be treated semiclassically and the quantum versions of the response 
and correlation functions will play important roles. 



III. QUANTUM FORMULATION 
A. The Fey nman- Vernon Approach 

In this section we describe the dynamics of the coupled oscillator from a quantum point 
of view. In order to do that we need, like in the classical case, a systematic way to eliminate 
the detailed information we don't need about the chaotic system. We will do that using the 
Feynman- Vernon approach Because of the non-linear chaotic system, we will not be 
able to perform exact calculations. Instead, we will resort to semiclassical approximations. 

We consider the quantum version of the full Hamiltonian, Eq.(^, and again we denote 
by g = (x, y) the pair of coordinates of the chaotic system. The density matrix operator can 
be written as 

p{T) = mT)){i;{T)\ = e-^^^/^|V^(0))(^(0)|e^^^/^ 

where ip{T) is the wave function of the whole system. In the position representation 

piz{T),qiT),z'iT),q\T)) = (z, g|V^(T)) (^(T) |/, g') 
= J dz(0)d/(0)dg(0)dg'(0)(2,g|e-*^^/^|2(0),g(0))(^(0),g(0)|^(0)) x 

(^(0)1^(0), g'(0))(z'(0),g'(0)|e*^^/^|z',g') 
= j dz(0)dz'(0)dg(0)dg'(0)ir(^(T),g(T),z(0),g(0))^(^(0),g(0))x 

K*{z'{T), q'{T), z'{0), g'(0))^*(z'(0), g'(0)), (26) 

where the propagators can be written in terms of Feynman path integrals as 

K(z(T),g(T),^(0),g(0)) = J Dz(t)Dg(t) exp 

with action 



^S[zit),qit)] 



S[z{t),q{t)] = r dt{Lo{z{t)) + LMt)) + Li{z{t),qm 
Jo 



(27) 



= So + S, + Sj. (28) 



Thus, 



p{z{T),q{T),z'{T),q'{T)) = 
= j d^(0)d^'(0)dg(0)dg'(0)D^(i)D^'(i)Dg(i)Dg'(i) 



exp 



{S[z{tlq{t)]-S[z'{tlcim 



p(;2(0),g(0),y(0),g'(0)), (29) 



where 

p(^(0), g(0), ^'(0), q'm = ^(^(0), g(0))r (^'(0), g'(0)) (30) 

is the initial state. As usual we use z{T) and q{T) just as labels for the positions z and g at 
T. z{T) and g(r) are not functions of T. 

We assume that the initial state can be written as 



p(;2(0), g(0), ^'(0), g'(0)) = Po{z{Q), y(0))p,(g(0), g'(0)) 
and define the reduced density matrix by 

Po{z{nz\T))^ j dq{T)piziT),q{T),z'{T),qiT)). 



(31) 



(32) 



We obtain 



Po{z{T),z'{T)) ^ J dz{0)dz'{0)Dz{t)Dz'{t) 

I J dq{0)dq'{0)dq{T)dq\T)Dq{t)Dq'{t)5{q{T) - q'{T)) 



exp 



{SMt)\ - So[q'it)] + S:[z{t), qit)] - S:[z'{t), q'{t) 



Pe(g(0),g'(0))^exp 
which can be written as 



h 



{sMt)]-So[z'm 



Po{z{0),z'{0)), 



(33) 



Po{z{T), z'{T)) = / d^(0)d^'(0) J(^(r), z\T), ^(0), ^'(0))p„(^(0), ^'(0)), (34) 



where 



J{z{T), z'{T), z{0), z'm = j Dzit)Dz'{t)J^[z{t), z'{t)] exp[-(5„[z(i)] - So[z'm] (35) 
is the 'superpropagator' and 

J^[z(t),z'(t)]^ J dq(0)dq'(0)dq(T)dq'(T)Dq(t)Dq'(t)5{q(T)-q'(T)) 



exp 



- {SMt)] - ScWrn + Sj[z{t), q{t)] - Si[z'{t),q'{t)]) 



Pe(g(0),g'(0)), (36) 



is the so called 'influence functional', which contains all the information about the chaotic 
system. 

Equation (|34j) is the equation of motion for the reduced density matrix. For = 1, 
J becomes the propagator of the isolated oscillator. Our goal is to get an approximate 
expression for J that includes the effects of the chaotic system. 

We take as the initial state for the chaotic system one of its energy eigenstate 0a('?) = 
(g|a). Thus, 

p,(g(O),g'(O)) = 0,(g(O))0:(g'(O)). (37) 

This is the quantum version of the classical microcanonical distribution we considered in 
section II. 

The difficulty in the calculation of J is that the chaotic Lagrangian Lc is not quadratic. 
Therefore, T has to be treated in a perturbative manner. Rewriting T as 



:F[z{t), z'{t)] = / dg(0) . . . Dg'(t)5(g(T) - g'(T)) exp 



(SMt)] - sAq'm 



exp 



&t{z{t)x{t) - z'{t)x'{t)) 



0a(g(O))0:(g'(O)) 



(38) 



we assume that 7 is small enough so that the exponential in the second line can be expanded 
to second order in its argument. These terms can be calculated by inserting complete sets 
of energy eigenstates of H^. The result, following |19i], is 



:F[z{t),z'{t)]^l-['^]x. 



dt[z{t) - z\t)] 



where 



Fa{t - S) = ^ exp[zu;fea(t - S)], CJfca 



Eh — En 



(39) 



(40) 



l{q)x(t)a{q), 



:F[z{t), z\t)] ^ 1 - \^W), z\t)] ^ exp 



(41) 



and El, are the eigen-energies of chaotic system [x is the coordinate oi Hc 'uiVi ). For the 
NS HJ—x) = HJx) and x„„ = 0. Thus, 



(42) 
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where 

,2 rT rt 



<l>[z,z'] = ^ dtj^ ds[z{t) - z'{t)][z{s)F:{t - s) - z'{s)F,{t - s)]. (43) 
With these approximations the superpropagator can be written as 

J{z{T),z'{T),z{0),z'{0))= [ D^(t)Dz'(t)e^(^«/W*)'-'Wl), (44) 



where we have defined the effective action 

S,j[z{t),z'{t)] = SMi)] - So[z\t)]+t^z{t),z'{t)]. (45) 

Since 5*6/ is quadratic in z and z' the path integral can be solved exactly by the stationary 
phase method. It is convenient to define the new variables r(t) = {z{t) + z'{t))/2 and 



y{t) = z{t) — z'{t) 2m ^^"i to separate Fa{t) = F^{t) + iF'^it) into real and imaginary parts. 



This allows us to write 

J(r(T),|/(T),r(0),y(0)) = j Dr(t)Dy(t) e^^l^^*)'^^*)] " ^'^['^^■^^1 (46) 

where 

T ( rt 



S[r{t),y{t)] ^ dt ^m[r{t)y{t) - uj',r{t)y{t)] + 2^'y{t) dsF^ {t - s)r{s) j , (47) 
and 

0[r(t),y(t)] = 7' Tdt fdsy{t)yis)Fiit-s), (48) 
Jo Jo 

are the real and imaginary parts of Sef- In Appendix A we show that 

J(r(r), y(T), r(0), 1/(0)) = G{T, 0) exp | ^S[r,, } exp { -^0[l/e, Ve] } , (49) 

where G{T,0) can be obtained by the normalization condition of reduced density matrix 
and re(t) and ye{t) are the extremum paths of S, which satisfy 

fe(t) + CUoVe(t) - / dsFj'it - S)reis) = 0, (50) 

fn Jo 

Ut) + ujMt) - ^ / dsFa"is - t)y,is) = 0. (51) 

Jo 

Therefore we need F^ to solve and and we also need F^ to calculate (j)[ye,ye\- 
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From Eq.dlO)) it follows that 



= (52) 



where x{t) is the Heisenberg representation of x. The real and imaginary parts of Fa are 

_ {a\{x{0),x{t)}\a) 

and 

„ (a|[:r(0),:r(t)]|a) 
^'^^'^ = 2^^^ ' (^^^ 

where {.} is the anticomutator and [.] is the comutator. Thus, F'^ and F^ are, respectively, 

the quantum analogs of the classical correlation and response functions of Section II. 



B. Semiclassical Expressions for Correlation Functions. 

In this section we obtain semiclassical formulas for F'^(t) and F'^{t). We write 

Kit) = ^Y.^h\{m.m}m = ^TrUm, (55) 

b 

where p = \a){a\ is the microcanonical distribution of the chaotic system and / = 
^x(0),x(t)\. To calculate F" we take / = —i[x{0),x{t)]. Using the Wigner-Weyl repre- 
sentation j25[ the trace can be written as 

Trlfm = J ^^fil^P-^t)W{q,p) (56) 

where 

POO 

f(q,p;t)= / dueiP-'^{q-u/2\m\q + u/2). (57) 



is the Weyl transformation (or symbol) of f{t) and W{q,p) is the Wigner function of p. For 
fit) = —i[x{0),x{t)] we have 



f{q,p;t) = due'P-'''/^ {q - u/2\{xx{t) - x{t)x)\q + u/2) 

J — oo 

POD 

= -I dMe*P-"/^(-u^)(g-M/2|x(t)|g + u/2) 



oo 

oo 



h^l I due'P-''/^{q-u/2\x{t)\q + u/2) 



dPx W-oo 

fOO 



due^P"/^ / dvVxK*{v,q-u/2;t)K{v,q + u/2;t)\- (58) 
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where u and v represent coordinates of the chaotic system and K{v,v';t) = {v\e~^^'''^^^\v') 
is the propagator. 

We now replace the propagators by their semiclassical expressions K and do the integrals 
by the stationary phase approximation. The stationary phase condition shows that the most 
important contributions come from the trajectories starting at q — u/2 (for K*) and q + u/2 
(for K) and arriving at v in the time t such that 

VMv, q - u/2) - V,Ri{v, q + u/2) = 0, (59) 

where Rk and Ri are Hamilton's principal functions coming from the phases in K and K*. 
Since VyRi{v, v') gives the final momentum, ()59|) imposes that the final momenta of the two 
trajectories must be equal. Since the final positions are also equal, the two trajectories must 
be identical. Thus, 

/oo 
dvv^k*{v,q-u/2;t)k{v,q + u/2;t) ^ x{q,p;t)6{u), (60) 
-oo 

where x{q,p;t) is the coordinate x of the stationary trajectory. 
Using (jnH in (jSHI) we find 

f{q,p;t) = h^x{q,p;t) = h{x{0),x{t)}. (61) 



smce 



J (()] (f]\ = ^^(0) ^^(^) _ ^3^(0) dx{t) ^ dx{t) 

^^>^^>^ dx{Q)dpM dpMdxiid) dpM ' 



For f{t) = {x{0),x{t)}, we find 



— oo 



f{q,p;t) = I duexp[-p-u]2q,{q-u/2\X{t)\q + u/2) 



^ 2q^x{q,p;t) = 2x{0)x{t). (63) 
The semiclassical limit of the Wigner function 

W{q,p;E) = ^j dte'^'^^ j duexp(j^p-u^ K{q + u/2,q-u/2;t). (64) 



was obtained by Berry !2g] and can be written as 

W{q,p; E) ^ 6{E - H{q,p)) + Wi{q,p- E). (65) 



13 



The first term is the classical micro-canonical distribution and the second, Wi, is given by 
classical periodic orbits corrections to the classical function. These periodic orbits have 
energy E = Ea, corresponding to the eigenstate \a) of the microcanonical quantum distri- 
bution. Using (j65|) we write 

1 



Tr[fit)p] 



dqdpf{q,p;t)6{E - H{q,p)) + / dqdpf{q,p;t)Wi{q,p;E) 
1 



[f + n ■ 



(66) 



When f{q,p] t) is replaced by the semiclassical expressions for the anticomutator and comu- 
tator, the first term of becomes, except for a normalization, the classical expressions 
for the response and correlation functions respectively. Following 2^ the second term of 
(IM)) becomes 



f\E,t)^J2A,cos{S,{E)/h + ^^) jdrf{q^{T),p^{r)-t). 



(67) 



where, Aj depends on the stability of j-periodic orbit, Sj{E) is its action, 7^- is the Maslov 
index and the integral is calculated over a period of the j-orbit. Analogous results can be 
found, for example, in |27|. 
Furthermore 



Trip) 



1 

(2^ 



j dqdp5{E - H{q,p)) + nq{E-h) 



n,{E) + n^{E-h)] 



(68) 



where the first term is the classical density of states and the second term is known as 
Gutzwiller's trace formula. Thus, 



Tr[f{t)p] _ f{E-t) + f\E;t) 



Tr[p] 



n,{E)+n,{E;h) 



(69) 



~ {fiq,p;'t))d. 

where {f{q,p;t)),i. = f%E,t)/n,{E). 



nq{E;h) 



+ 



f\E;t) 
nc{E) 



nq{E;h) 
na(E) 
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Finally, we can calculate the semiclassical expressions for F^'(t) and F'^{t). From dHU), 
dnSl) and (|69|) 



Kit) 



(X(0)x(t))e/. 



_ nq{E;h) 
nc{E) 



(70) 



and 



(0(t)) ({x(0),x(t)}),. 



2/1 



_ ng{E; h) 
nc{E) 



+ 



2ra, 



^ A, cos + 7,) dr{x,(r), x,(r + t)}. 



(71) 



Both these semiclassical expressions are given by their classical counterparts multiplied 
by a correction to their amplitudes, given by ng{E; h)/nc{E), plus a correction from periodic 
orbits. The temporal dependence of the first is given solely by the classical dynamics and 
it decays exponentially. The second term, however, is a sum of oscillating functions and 
carries the temporal dependence characteristic of the chaotic system. As a final remark we 
note in Q| these functions were calculated using random matrix theory. 



C. The Superpropagator 

The semiclassical expressions for F'^{t) and F'^{t) allows us to solve the equations of 
motion and and to calculate the superpropagator, Eq. ipUl) . In order to have 
explicit formula, we shall consider only the zero order approximation F^(t) ^ (f)xx{t)/2. 
Substituting it in (jSUj) and and integrating by parts we get 

fe{t) + Qlre{t) + ^ [ dsF{t-s)r4s) = -—r,{0)F{t), (72) 

i/e(t) + Xo2/e(^)-- / dsFit-s)Us) = ^2/e(0)F(t), (73) 
where = uJq — 7^F(0)/m, Xo = + 7^-^(0) /m and 

F{t) = j ds<PUt-s), (74) 

as in the classical case (eq.(|^). 
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In Appendix B we solve these equations by the method of Laplace transforms. We find 
that the solutions involve two very different time scales. The shortest time scale is relevant 
only for times much smaller than the period of the decoupled oscillator. For times 

of the order of 1/ujq these terms can be discarded as transients. Here we shall adopt this 
approximation and keep only the terms that are significant for times of 1/ujq. In this case 
we show that ()72j) and ()73p can be written as 

h{t) + 2Are(t) + ^lre{t) = 0, V?^:^ K\ (75) 
Veit) - 2AUt) + xlVeit) = 0, xo > A', (76) 



where 



A = — lim f dsF(t-s). (77) 
Approximating ~ Xo ~ we find 

. I Sin Mr- 01 ^( ^ e-i^^.(T) 

sm [ojqT) sm [ujqT) 

1^ sm [ijJqT) sm [uJqT) J 

Within this approximation we can calculate the real part S of the effective action, Eg. ()47j] . 
We obtain 

S[r^, Ve] = [rmolK{T) - mA]r{T)y{T) + [mujlK{T) + mA]r{0)y{0) 

-mwlL{T)riQ)y{T) - mu;',N{T)r{T)y{0), (79) 



where, 



= (80) 

ujq sm[uJoI ) 



-AT 



L{T) = — --, (81) 

ujq sm[uJoI ) 

1 e^^ 

N{T) = --^r^y (82) 

ujq sm[ujQl ) 



For the imaginary part of the effective action, Eq. pHj) . we get 



<P[ye,ye\=l' [ dty,{t) [ dsF;^{t-s)ye{s). (83) 
Jo Jo 
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Using again only the zero order term of the semiclassical expression dZOl) for and the fact 
that the classical correlation function decays exponentially, as in (fTTjl . we can approximate, 
for t ~ l/uJo 

B' 



I 



t rt 

dsF'^{t-s)y,{s)^y,{t)\im / ds Fi{t - s) = y^it)- 



t—*oo 



which implies that 



2 Ul rT 



dtylit) 



[A{T)y\T) + B{T)y{T)y{Q) + C{T)y' 



(84) 



(85) 



where 



A{T) 
B{T) 
C{T) 



+ (A^ + uol) - cos {2uoT) - Aujq sin (2woT)] 



4A(cu2 + A2) sin^ (cuoT) 
[—uJq sinh (AT) cos {ujqT) + Aujq cosh (AT) sin (cuqT)] 
A(w2 + A2) sin^ (uJoT) ' 
^^2g2AT _ ^ ^2^) _^ ^2 (2cjoT) - Acuq sin (2woT)] 



Finally, putting everything together we get 

J(r(T), 2/(T), r(0), 1/(0)) = G(T, 0) exp | ^S[r,, y,] | exp |-^0[l/e, 1/e 
= G(T,0)exp|^i^2(T)r(T)y(T)|exp|-ii(T)y2(T)| x 



exp 



exp 



n 



K^{T)r{0)y{0) - L{T)r{0)y{T) - N{T)r{T)y{0) 



X 



^^B{T)y{T)y{0) - \c{T)y\0) 



where 



i^i 2(T) = muolK{T) ± mA L(T) = mcu^L(T) iV(T) = mulN{T) 



Am = i AiT) BiT) 

n 



7^5' 

5(T) C(T) = -— C(T). 

n 



(86) 
(87) 
(88) 



(89) 



(90) 



We finish this section with a comment about the physical situation described by these 
calculations. Since we have considered only the first terms of the semiclassical expressions 
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for F'^ and F'^ our results are valid only in the Ehrenfest time scale of the chaotic system 
given by [3| 

J- .„(!). (91) 

where is the Lyapunov exponent and Sc is a typical action of the chaotic system, for 
example the action of the shortest periodic orbit. Approximating Al by a, we see that, in 
order to observe effects for t ~ I/ujqi we must have 

i,.i.i|„(|)^|.expfi^). (92) 

For NS, a/cuo ~ 8 and Sc ~ 10 which means that h must be smaller than 10~^ for our results 
to be valid. 



IV. APPLICATIONS 



The superpropagator allows us to study the time evolution of the oscillator under the 
influence of the chaotic system. The reduced density matrix satisfies 

Po{riT),y{T)) = J dr(0)di/(0) J(r(T), y(T), r(0), i/(0)) p,(r(0), y{0)). (93) 

In the following we calculate explicitly the propagation of two different oscillator's states. 
These two applications are similar to the ones presented by Caldeira and Leggett to study 
dissipation and decoherence. 



A. Propagation of a Gaussian State. 

For a Gaussian state 

mO)) = j^^,e'nMoy^e-^Hoy^'^' (94) 

the density matrix po(-2(0), -2'(0)) = ip*{z'{0))ilj{z{0)) can be written in terms of r = {z+z')/2 
and y = z — z' as 
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Substituting (jH^ and in and performing the integrals, we get 



p,(r(T),y(T)) = G(T,0) 
X exp 

X exp 



27rr 



1/2 



2[2aC'i(T)+(T2i<-2(T)] 



r(T)- 



Af(r) 



A(T) 

h 



2^^ 



{a2Ki{T)L(T)-hB(T))^ 



xe^p\MT)riT)yiT)-^ 



(cT^KiiT)L{T)-hB{T)) 



iV(T) fr(T) 



Ar(T) 



(96) 



where, 



CiiT) = C{T) + 



n 



8cx 



2 ■ 



(97) 



For y = z — z' = 0, Po becomes the probabihty density. After normahzing we obtain 

1/2 

Po{r{T),0) 



X exp 



27r[2hCi{T) + a^KfiT)] 



r{T) 



P 



Eq. 



2[2hCi{T) + a^Kf{T)] \ ' ' N{T) 
represents a Gaussian packet whose center follows the trajectory 



(98) 



r(T) 



V 



V 



e sin(ci;o^)- 



(99) 



N{T) mujo 

The dissipative effect due to the interaction with the chaotic system is explicit. The same 
behavior was obtained by Caldeira and Leggett ^ using a thermal bath with many degrees 
of freedom. Eq.()99p represents the trajectory of a weakly damped harmonic oscillator. The 
critical and strongly damped cases cannot be described by this formalism because of the 
weak coupling regime adopted. 

The width of the evolved packet is given by 

a^kf{T) + 2hCi{T) 



a\T) 



m{T) 



(100) 



After some algebra we can show that 
AT) 



a 



1 + e2)e-2A^ 



l + e2 



ril - e-2Ar 



;i + 2e sin (luoT) cos (cuqT) + 2e^ sin^ (cuoT))] 1 

^ TT7^ J' 



101) 
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where 

..A r = ^. (102) 

The expression above for F comes from the following considerations: from Ci{T) it follows 
that 



Using Eqs.((771) and (|2S1) and the relation [24 1 



(103) 



(P..(0)x(t))e = -^(x(0)a;(t))e (104) 



we obtain 



7^5' 

mE,(0)' ^ ^ 

which leads directly to p02|) . Notice that, due to (j^ . F ^ 1 is the only possibility. 

Fig. 3 shows that cr^(T) for F = 1, 0.5 and 2.0. These curves can be well fitted by the 
simpler expression 

a\T) = a2[e-2^^ + F(l - e'^^^)], (106) 

which, for t ~ l/ciJo, can be written as 

a^{T)=a^[l + {T -l)2KT]. (107) 

We see that (7^(F — 1)A plays the role of a diffusion constant. Fig. 3 also shows that F controls 
the increase or decrease of a^{T). In the present case, cr^(T) can only increase because of 
the constraint F ^ 1. In the Caldeira-Leggett model, on the other hand, the width can also 
decrease if the the temperature is very low. 



B. Superposition of Two Gaussian States. 



We now consider an initial state consisting of two Gaussian wave-packets, one at the 
origin and one centered at z{0) = q^: 



7^(^(0)) = ArV2[^^(;,(o)) +^2(^(0))] 

z\m , r (^(o)-go)' 



exp 



4a2 



+ exp 



4a2 



(108) 
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The density matrix is given by 

Po{z{^),z'm= N [pn(^(0),z'(0))+p22(^(0),/(0)) 
+ pi2(^(0),/(0))+p2i(^(0),^'(0))]. 



(109) 



with pij = The time evolution of po can again be calculated with Eq.(|^. The 

result, for y = z — z' = is 

1/2 



PiMT),0) 



P22(r(T),0) 




2[1 + MT)] \nf{T) 
pi2(r(T),0) + P2i(r(T),0) 



exp 



exp 



fiT) 
N^jT) 



'r{T)-Q{T)f 



(110) 
(111) 



N' 



;(T)\ 



1/2 



2[1 + MT)] \T:f{T)l 



exp 



8a' 



9iT) 



exp 



X exp 



X 2 cos < ■ 



fiT) 
hN\T) 



[4a^f{T)K^{T) 



r{T)-Q{T)Y 

[(r(T) - Q{T)f - r\T)\ 



where 



/(T) = 2[2hC,{T)+a'Ki{T)] 



h{T) = exp 



[l + ^7(r)] 



Q{T) = ^— — go, 



A^(T) 



2hC{T) 



2nCi{T) + a^Kl{T) 
The interference term can also be rewritten as 

pi2(r(r),0) + P21 (r(T), 0) = 2 COS [a(T)((r(T)-g(T))2-r2(T))] 



(112) 

(113) 
(114) 

(115) 

(116) 



X pf (r(T),0)p2'^(r(T),0)exp 



1/2, 



9{T) 



(117) 



Eq. ()117|) shows that the interference is attenuated by exp [— (q n/8( 7^)qfT)]. Eq. piTj) is 



very similar to the expression obtained by Caldeira and Leggett 
temperature dependence in g{T)^ which can be written as 

Th{T) 



14| . although there is no 



9iT) 



:i + e2)+r6(T)' 



(118) 
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with 



b{T) = e^^^ - 1 - 2e sin (cuoT) cos (cuqT) - 2e^ sin^ (cjqT) 



;il9) 



and e = A/ujq. We note that the asymptotic hmits 

^(T = 0) = and g{T ^ oo) ^ 1 (120) 

are the same as those in the Caldeira-Leggett model. 

Fig.4 shows g{T) for F = 10. In the regime F ^ 1, we can approximate g{T) by 

2FAT 

9iT) = (121) 



1 + 2FAT' 

This simphfied expression helps to estimate of the decoherence time. For example, with 
(11211) . we can estimate the time T' such that 



exp 



10" 



;i22) 



Defining n = Q'o/Scr^ (the number of quanta Huq of the wave packet centered at go); "we get 



n-ln(lO) 
31n(10) 



2FAr' = 2nFAr' = 1 ^ T' 



1 



(123) 



2nFA 

Since we are interested in the situation where n ^ 1 and F ^ 1, we find that the decoherence 
time is much smaller than the time scale where dissipation takes place, i.e., T' ^ 1/A. 



V. DISCUSSION AND CONCLUSIONS 



We have made two important assumptions in our calculation of the superpropagator. 
The first of these assumptions, the weak coupling regime, was important to reduce the path 
integral to a quadratic form in the oscillator variables. The second assumption was the 
semiclassical regime of the chaotic system. This was essential to establish the connection 
between the coupling in the influence functional and the classical correlation and response 
functions that enter in the classical description of the system. The use of these classical 
functions make the importance of the chaotic dynamics explicit and show that the time 
scales obtained classically are important ingredients to describe dissipation. In particular, 
the exponential decay of correlations happens in a time scale much shorter than the natural 
period of the oscillator. The time of correlation loss plays the role of the microscopic time 
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scale in the Brownian motion, which is much shorter than the macroscopic one [29[. More- 
over, the exponential decay of the classical correlations is what makes dissipation possible in 
the present treatment. The corrections due to periodic orbits have not been explored here 
and the importance of their contribution to dissipation and decoherence is not clear at this 
point. 

The effective dynamics we obtained, expressed in (jH^ . is analogous to the Caldeira- 
Leggett theory in the limit of high temperatures and weak damping jisl . Q| . For example, 
the diffusion constant in ()107|) can be written, for F ^ 1, as 



EM 

cr i ii = 

which should be compared with 



a'TA = ^^A, (124) 



D = (125) 

mujQ 

for the Brownian motion. Therefore, -E'c(O) plays the role of fc^T. From Fig. 4, F seems to 
play the role of ksT /huok since it controls the behavior of ct^(T). However, despite this close 
analogy between the two models, our results are valid only for short times since they are 
limited by Ehrenfest time and perturbation theory. 

In summary, we have shown, using Feynman- Vernon approach, that a chaotic system with 
two degrees of freedom can induce dissipation and decoherence in a simple quantum system 
when weakly coupled to it. The formalism we have chosen allows us a close analogy with 
the many body formulation of the Caldeira-Leggett model. The most important quantities 
in the formalism, the correlation and response functions, are obtained directly from the 
dynamics, and not from phenomenological assumptions as in the Caldeira-Legget model. In 
our approach we have used simple classical approximations and discarded all periodic orbits 
corrections. The effects of these corrections are certainly worth studying. 



APPENDIX A: THE STATIONARY PHASE APPROXIMATION 

In this appendix we solve the path integral Eq. ()46p by the stationary phase approxima- 
tion. Let {rf>{t),ye{t)) be the stationary path and 

r{t) = r,{t) + 5r{t) = r,{t) + eif(t) (Al) 
y{t) = ye{t) + 6y{t) = ye{t) + e^m , (A2) 
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be a neighboring path with r(T) = r(0) = and y(T) = y{0) = 0. 
The stationary path is obtained from the condition 

AS = S[re{t) + e,r{t),yeit) + e2y{t)] - S[re{t),ye{t)] 

dAS I ^ dAS 



+^2^ - U 



We find 

dAS 



dei 

and 

dAS 



de-2 



^ dtf (t) |m[i/e(t) + - 27^ ^ dsF/(s - t)l/e(s)| , 



dty{t) \m[f,{t) + cUoVe(t)] - 272 f dsFj\t - s)r,{s) 
I Jo 



where we used f(T) = f(0) = 0, y{T) = y{0) = and dt ds = 2 dt /J ds. 
Therefore, the equations of motion for the stationary path are given by 

„2 



2'Y r 

Ut) + i^heit) - — / dsFa"{s - t)ye{s) = 

m Jo 



and 

„2 rt 



hit) + CUoVe(t) - ^ / dsFa"{t - S)re(s) = 0. 

m Jo 

Expanding 0, Eq HHl around the stationary path we find 



[re(t) + eif(t),ye(t) + e2y(t)] 



= ^7^/ dt [ ds[y,{t) + e2ymye{s) + e2y{s)Wa{t-s) 

^ Jo Jo 

= ^7^ [ dt [ dsyeit)yeis)Fa'it -s) + 

^ Jo Jo 

\l^^2 I dt I ds[y{t)ye{t)+ye{s)y{s)]F'a{t-s) + 

^ Jo Jo 



^2ifi I / dsy(t)y{s)F^{t- s) = <i)[ye,ye] + 2'^[y,ye] + '{'[y,y\■ 



'o JO 
Therefore, from P6|l . we have 



J(r(T),y(T),r(0),y(0)) = 

/•o 
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We are now going to show that ()A9|) is a function of the initial and final times only, 
which is not obvious because of the functional dependence on ye{t). In order to do this we 
discretize the paths and re-write ()A9|) in the form 



exp <j ^S[6r, Sy] 



N 

exp 



N N 



3=1 k=l 



(AlO) 



where Sr^ = 6r{tj), F^^^_^^ = {tj - 4), 

^ I i=l k=l ) 



and 



exp !^-lv[Sy, ^ exp E^^'e' E ^%-^-^^<^o-., | • (^^2) 



Grouping the exponents we obtain 



exp!^^S[6r,6y] - ^v[Sy,Sy] - ^^[Sy,yS^ ^ exp !^-'-U^MU - A^U^, (A13) 



with 



U'^ = {6n...6rN6y,...6yN) M=| °^ |, (AM) 

p r 

and where p and r are NxN matrices and = (0 a) and a are A^- dimensional 
vectors. To solve the path integral we need to integrate this exponent over dU 
d6ri . . . d5rjv d6yi . . .Sy^. The result is 

1 I 1 



^exp ^-\a^M-'A 



(detM)i/2-^ I - - ■ ^^'^^^ 

Because M has a zero upper left block, its inverse has a zero lower right block and, 
therefore, A'^M~^A = 0. Since all the dependence on the initial and final positions is 
contained in A, ()A9|) is indeed a function only of the initial and final times. Therefore we 
may write the superpropagator as 

J(r(r), |/(T), r(0), 1/(0)) = G{T, 0) exp | '-S[r,, } exp | -^0[l/e, ye] } , (A16) 
and G{T, 0) can be calculated by imposing the normalization of the reduced density operator. 
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APPENDIX B: SOLUTION OF THE EQUATIONS OF MOTION. 



Taking the Laplace transform of (j^Uj) . we get (with F'^{t) ~ (f)xxi't)/2) 



m 



re(s) = sr(0) +r(0). 



where f{s) = C{fit)} is the Laplace transform of f{t). Using 

2 



-(Px(0)x(t))e = Ae-°l*l sin (ut), 



BTI) becomes 



rAs 



s[{s + + uj^]r{0) + [{s + af + u;^]r{0) 



(Bl) 



(B2) 



(B3) 



The Heaviside's theorem establishes that if P{s) and Q{s) are polynomials such that the 
order of P{s) is smaller than the order of Q{s), then 



-1 



Pis) 

Q{s) 



y- 



PisA 



(B4) 



where Si are the roots of Q{s) = and Q'{s) is the s-derivative of Q(s). Therefore we need 
the roots of 



X 



a 



[x + ir + {- 

.a 



7^ Au 



m 



0, 



where x = s/a and Qq ~ Uq. From Section II we have 



— V ~ 1.6 X 10"^ 
a 



a J 
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7^ Auj o 
I _ ~ 3 X 10-^ 



m 



(B5) 



(B6) 



and the roots of (IB5I1 are 



xi = -1.00 -iS.OO 
^3 = -4 X 10"^ - z0.12 

Multiplying these roots by a, we get 

Si ~ —a — iuj 
S3 -A - iLJo 



X2 = -1.00 + i5.00 

a;4 = -4 X 10-^ + 20.12, 



(B7) 



S2 ~ —a + iuj 
Si ~ —A + iujQ. 



The same procedure is applied to (j^J)- The Laplace transform of (jSlj) is written as 

s[{s + af + a;%(0) + [{s + af + uj^]m 



ye[s) 



{{s^ + Ql)[{s + af + uj^] + ^^Au;} 



(B8) 



(B9) 
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and the roots are 



Si ~ —a — luj S2 ~ —a + luj 

S3 ^ A — iu}Q S4 ^ A + iuQ. (BIO) 

Since we are interested on time scales such that t ~ I/cuq, Si and S2 are transient solutions 
and only S3 and S4 are important. Therefore, turning to the equations (f7^ and ()73p and 
considering times on the scale t ~ I/ojq, we see that those equations can be rewritten 
approximately as 

r,{t) + 2Ar,{t) + nlr,{t) = Q, (Bll) 
i/e(t)-2Ay,(t) + Xoye(t) = 0, (B12) 

where terms proportional to F{t) were disregarded (since they go to zero for t ~ I/cuq) and 
the convolutions terms were approximated in the following way 

/ dsF(t- s)re(s) ~ re(t) lim / dsF{t-s). (B13) 
Jo io 

Thus, A is given by 

A = — lim f dsF{t - s). (B14) 
2m t^oo Jq 

Indeed, applying the Laplace transform in (jBlljl and ()B12jl . we get the roots 

si = —A — iuo S2 = —A + iuo, (B15) 

for re{t) and 

si = A — iujQ S2 = A + iuo, (B16) 



for ye{t) since ^l,Xo > A and fig ^ Xo ~ ^o- Comparing ()B15jl and ()B16j) with (|B8|l 
and (jBlOj) . we conclude that the equations (jBllj) and (jB12j) give a good description of the 
behavior given by (ff^ and (ffH|) for t ~ l/c^o- 
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FIG. 1: Correlation functions for the NS for = 0.38: (a) {p^{0)x{t))e; (b) {x{0)x{t))e. The full 
line shows the numerical results and the dashed line shows the fitting. The averages were computed 
using 35000 initial conditions. 
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FIG. 2: Average oscillator energy at short times with the NS as chaotic system. Tq = I/loq. 
The dashed line shows {Eo(t)) and the doted line shows {Eor{t)), both obtained numerically. The 
full line shows Eq.^ without f{t). (a) EoiO)/Ec{0) = 1.0, (b) EoiO)/Ec{0) = 0.25 and (c) 
Eo{0) / Ec{0) = 0.1. The oscillator's parameters, coupling constant and number of initial conditions 
are the same as in Fig. 2. 



29 



1.006 



1.004 



1.002 



1.000 



0.998 



0.996 





1 T/T„ 



12 16 




FIG. 3: Squared width of wave packet cj^(T)/o"^ as given by Ea. ()lUl() with Tq = I/loq. (a) T = 1.0, 
(b) r = 0.5 and (c) T = 2.0. 
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FIG. 4: The full line shows ^(r) as given by Ea. lfTTH|) for T = 10.0. Tq = I/loq. The dotted line 
shows g{T) as given by Ea. (|121() for F = 10.0. 
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